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Abstract 

The RG expansions for renormalized coupling constants and of the 
3D n-vector model are calculated in the 4-loop and 3-loop approximations 
respectively. Resummation of the RG series for g^ by the Pade-Borel-Leroy 
technique results in the estimates for its universal critical values gUE) which 
are arqued to deviate from the exact numbers by less than 0.3 %. The RG 
expansion for g% demonstrates stronger divergence being much less suitable 
for getting proper numerical estimates. 

Recently, the free energy (effective action) and, in particular, higher-order renor¬ 
malized coupling constants g 2 k for the basic models of phase transitions became the 
target of intensive theoretical studies (see, e. g. ii and references therein). These 
constants enter the scaling equation of state and thus play a key role at criticality. 
Along with critical exponents, they are universal, i.e. possess, under T ^ Tc, nu¬ 
merical values which depend only on the space dimensionality and the symmetry 
of the order parameter. Calculation of the universal critical values of gs, etc. 
for the 3D Ising model by different methods showed that the field- theoretical RG 
approach in hxed dimensions yields the most accurate numerical estimates. It is a 
consequence of a rapid convergence of the iteration schemes originating from RG 
expansions iii- It is natural, therefore, to use the held theory for calculation 
of renormalized higher-order coupling constants for the 3D n-vector model. In the 
report, the 3D RG expansion for the renormalized coupling constants g^ and gs will 
be found and the universal critical values of g^ will be estimated. 
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The 3D 0(n)-symmetric model is described at criticality by Euclidean scalar 
field theory with the Hamiltonian 
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(TX 




( 1 ) 


where rug is proportional to T — TjP\ being the phase transition temperature in 
the absence of the order parameter fluctuations. The fluctuations give rise to many- 
point correlations < Lp{xi)Lp{x 2 )---^{x 2 k) > and, correspondingly, to higher-order 
terms in the expansion of the free energy in powers of the magnetization M: 

OO 

F(M, m) = F( 0 , m) + Y, , ( 2 ) 

fc=i 

where m is a renormalized mass, r; is a Fisher exponent, and g 2 k are dimensionless 
coupling constants. Let, as usually, g 2 = 1/2. Then g 4 , g^, fi'sv will acquire, under 
T —> Tc, the universal values. 

The asymptotic critical values of 5 ^ 4 , (74 (n), determining critical exponents and 
other universal quantitites, have been found from the 6 -loop expansion for RG (3- 
function i.S.1.0; they are known with rather high accuracy. On the contrary, 
the information about the universal numbers gl, g^, etc. for 0(?7,)-symmetric model 
is very poor today. To estimate them, we calculate corresponding RG series and 
perform their resummation by means of the Pade-Borel-Leroy technique. The RG 
series for g^ and gs are obtained from conventional Feynman graph expansions for the 
6 -point and 8 -point vertices in terms of the bare coupling constant A. In its turn, A 
is expressed perturbatively via the renormalized coupling constant g^. Substituting 
then the series for A into the ’’bare” expansions, we obtain the RG expansions for 
ge and gs- 

As was shown |^, the 1-, 2 -, 3-, and 4-loop contributions to g^ are formed by 
1, 3, 16, and 94 one-particle irreducible Feynman graphs, respectively. In this case, 
the calculations just described give: 
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- (-0.00000949 + 0.00783129 + 0.34565683 n + 2.14825455)c/^ . 


ga = -gl 

71 


In the case of gs, the 1-, 2-, and 3-loop contributions are given by 1, 5, and 36 
graphs respectively p. Gorresponding RG expansion is found to be: 
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+ (0.00943497 + 0.60941312 n + 7.15615323) 
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Being a field-theoretical perturbative expansions these series are divergent (asymp¬ 
totic). To get reasonable numerical estimates for Qq and some procedure making 
them convergent should be applied. The Borel-Leroy transformation 


fix) = £ CiX* = / t^e ^F{xt)dt, F{y) = ^ 


i=0 
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can play a role of such a procedure. Since the RG series considered turns out to be 
alternating the analytical continuation of the Borel-Leroy transform may be then 
performed by using Bade approximants [L/M]. 

For qq we have the 4-loop RG expansion and can construct, in principle, three 
different Fade approximants: [2/1], [1/2], and [0/3]. To obtain proper approxima¬ 
tion schemes, however, only diagonal [L/L] and near-diagonal Bade approximants 
should be employed. That’s why further we limit ourselves with approximants [2/1] 
and [1/2]. Moreover, the diagonal Fade approximant [1/1] will be also dealt with 
although this corresponds to the usage of the lower-order, 3-loop approximation. 

The algorithm of estimating we use here is as follows. Since the Taylor ex¬ 
pansion for the effective action contains as coefficients the ratios i? 2 fc = 92 k/ 
we work with the RG series for Rq. It is resummed in three different ways based 
on the Fade approximants just mentioned. The Borel-Leroy integral is evaluated 
as a function of the parameter b under g^ = gl{n). For the hxed point coordinate 
glin) the values extracted from the six-loop RG expansion are adopted [jl|, ||. The 
optimal value of b providing the fastest convergence of the iteration scheme is then 
determined. It is deduced from the condition that the Fade approximants employed 
should give, for b = bopt, the values of Rg which are as close as possible to each 
other. Finally, the average over three estimates for Rg is found and claimed to be a 
numerical value of this universal ratio. 

The results of our calculations are presented in Table 1. It contains numerical 
estimates for g^ resulting from the 4-loop RG expansion (column 3) and their analogs 
given by the Fade-Borel resummed 3-loop RG series |]^ (column 4). As is seen, 
with increasing n the difference between the 4-loop and 3-loop estimates rapidly 
diminishes: being small (0.9 %) even for n = 1, it becomes negligible at n = 10 
and practically disappears for n > 14. Such a behaviour is quite natural since with 
increasing n the approximating properties of RG series for g^ are improving |l| . 

How close to the exact values of g^ may the numbers in column 3 be? To clear up 
this point, let us compare our 4-loop estimate for Rg at n = 1 with those obtained 
recently by an analysis of the 5-loop scaling equation of state for the 3D Ising model 
0, R. Guida and J. Zinn-Justin have obtained Rg = 1.644 and, taking into 
account some additional information, Rg = 1.643, while our estimate is Rg = 1.648. 
Keeping in mind that the exact value of Rg should lie between the 4-loop and 5- 
loop estimates (the RG series is alternating), our estimate can differ from the exact 
number by no more than 0.3 %. Since for n > 1 the RG expansion (^) should provide 
better numerical estimates than in the Ising case, this value (0.3 %) represents an 


3 



upper bound for the deviation of the numbers in column 3 of Table 1 from their 
exact counterparts. 

It is interesting to compare our estimates with those obtained by other methods. 
Since 1994, the universal values of the sextic coupling constant for the 3D 0{n)- 
symmetric model were estimated by solving the exact RG equations |^, by lattice 
calculations , and by a constrained analysis of the e-expansion ; corresponding 
results are collected in columns 5, 6, and 7 of Table 1 respectively. As is seen, they 
are, in general, in accord with ours. For large n, our estimates are consistent also 
with those given by the 1/n-expansion which are presented in column 8. 

The RG expansion for the octic coupling constant turns out to be much worse 
than the series from the point of view of their summability. Indeed, the series 

diverges considerably stronger and is one term shorter than that for qq. It 
implies that the only Fade approximant - [1/1] - may be really used in a course of 
the resummation of this series. In the Ising case n = 1, such a simple Pade-Borel 
procedure, when applied to the 3-loop RG expansion for gs, was found to lead to 
rather crude numerical estimates [^. As our analysis shows, with increasing n the 
situation becomes better but, nevertheless, the RG estimates for gl{n) remain much 
less accurate than those obtained for the sextic coupling constant. Gorresponding 
numerical results will be published elsewhere. 
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Table 1: Our estimates of universal critical values of the renormalized sextic coupling 
constant for the 3D n-vector model (column 3). The fixed point coordinates g* are 
taken from [5] (1 < n < 3) and [1] (4 < n < 40). The Qq estimates extracted 
from the Pade-Borel resummed 3-loop RG expansion (column 4), from the exact 
RG equations (column 5), obtained by the lattice calculations (column 6), resulting 
from a constrained analysis of the e-expansions (column 7), and given by the 1/n- 
expansion (column 8) are presented for comparison. 


n 

g* 

g*6 


»«• M 


96*0 

9*6 (IM) 


2 

3 

4 

5 

6 

7 

8 

1 

1.415 

1.608 

1.622 

1.52 

1.92(24) 

1.609(9) 


2 

1.406 

1.228 

1.236 

1.14 

1.27(25) 

1.21(7) 


3 

1.392 

0.951 

0.956 

0.88 

0.93(20) 

0.931(46) 


4 

1.3745 

0.747 

0.751 

0.68 

0.62(15) 

0.725(29) 

1.6449 

5 

1.3565 

0.596 

0.599 




1.0528 

6 

1.3385 

0.483 

0.485 




0.7311 

7 

1.321 

0.396 

0.398 




0.5371 

8 

1.3045 

0.329 

0.331 



0.319(4) 

0.4112 

9 

1.289 

0.277 

0.278 




0.3249 

10 

1.2745 

0.235 

0.236 




0.2632 

12 

1.2487 

0.174 

0.175 




0.1828 

14 

1.2266 

0.134 

0.134 




0.1343 

16 

1.2077 

0.105 

0.105 



0.1032(4) 

0.1028 

18 

1.1914 

0.0845 

0.0847 




0.0812 

20 

1.1773 

0.0693 

0.0694 




0.0658 

24 

1.1542 

0.0487 

0.0488 




0.0457 

28 

1.1361 

0.0360 

0.0361 




0.0336 

32 

1.1218 

0.0276 

0.0276 



0.0275(1) 

0.0257 

36 

1.1099 

0.0218 

0.0218 




0.0203 

40 

1.1003 

0.0176 

0.0176 




0.0164 
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